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1. Introduction 

One of the most appealing possible extensions 
of the Standard Model (SM) is supersymmetry 
(SUSY) [T]. It would stabilize the electroweak 
mass hierarchy and facilitate grand unification, it 
predicts a relatively light Higgs boson that would 
be consistent with the indications from precision 
electroweak data, it offers a possible explanation 
of the apparent discrepancy between the experi- 
mental measurement of the anomalous magnetic 
moment of the muon, (3 — 2)^, and the theoretical 
value calculated within the SM, and the lightest 
supersymmetric particle (LSP) is a plausible can- 
didate for astrophysical dark matter. 

Wc have published results from frequentist 
analyses of the minimal supersymmetric exten- 
sion of the Standard Model (MSSM), using like- 
lihood functions to take into account the experi- 
mental, phenomenological and astrophysical con- 
straints on SUSY. The unconstrained MSSM con- 
tains too many parameters for a full exploration 
of its parameter space to be possible using present 
data, even including the current LHC data set 
of ~ 35/pb [2]. Therefore, we have focused on 
making estimates within simplified versions of the 
MSSM, specifically the constrained MSSM (the 
CMSSM) [H] in which soft SUSY-brcaking mass 
parameters are assumed to be universal at the 
GUT scale, in the simplest generalization of this 
model in which the universality is relaxed to allow 
non-universal Higgs masses (the NUHMl) [SHI], 
in a very constrained model in which the supple- 
mentary relation Aq = Bq + mg is imposed on 
trilinear and bilinear soft SUSY-breaking masses 
in the CMSSM (the VCMSSM) [HE], and in min- 
imal supcrgravity (mSUGRA) in which, in ad- 
dition, the gravitino mass 7713/2 is set equal to 
the common soft SUSY-brcaking scalar mass mg 
before renormalization [8l[9]. In each case, we 
assume that the LSP is the lightest ncutralino 
x5- More details on the model definitions can be 
found in [8]. 

In a series of papers [HIBHS] we have presented 
predictions for Higgs and sparticle masses as well 
as for BR(i?s — >■ fi^ fi^) and the spin-independent 

^We recall that our convention [2l|8] for the sign of Aq is 
opposite to that of Sof tSUSY. 



dark matter scattering cross section, , and also 
for nit and M^r [TU]. Most recently [^ we have 
included in global analyses the results of an ini- 
tial CMS search in multijet + channels (CMS 
ax) [H] and an ATLAS search in lepton -I- mul- 
tijet + channels (ATLAS IL) [12] H Incor- 
porating these new results led to upward shifts 
in the lower bounds on the gluino mass, m^, 
by ^ 100 GeV in the models considered. Other 
masses connected to , such as that of the light- 
est neutralino, m^o (which we assume to provide 
the astrophysical cold dark matter (CDM) [H] ) 
also moved upward by corresponding amounts. 
This in turn led to somewhat lower expectations 
for the spin-independent dark matter scattering 
cross section cTp^ in the models considered [2]. 

Subsequent to our analysis [2] of the impli- 
cations of these initial LHC searches for SUSY, 
LHC experiments have provided several new con- 
straints on SUSY using an integrated luminosity 
of - 35/pb of data at 7 TeV. ATLAS has pub- 
lished the results of a search in multijet -I- I^t 
channels (ATLAS OL) [15] that has greater sensi- 
tivity in some regions to the types of gluino and 
squark pair-production events expected in the su- 
persymmetric models discussed here than did the 
earlier ATLAS IL search [12], and has also re- 
leased results obtained by combining the one- and 
zero- lepton searches [16] . CMS has announced re- 
sults from two other searches in multijet + 
channels that improve the CMS ax sensitivity 
also to gluino and squark production in the mod- 
els discussed here. Both ATLAS and CMS have 
also published the results of searches for jets -|- 
events with b tags [17], and for multilepton 
+ jets events [T^. In addition, CMS and 

ATLAS have published new upper limits on the 
production of the heavier neutral MSSM Higgs 
bosons H, A [l9l[20] , and LHCb has recently pro- 
vided a new upper limit on BR(i?s ^~^fi~) [21], 
of comparable sensitivity to previous results from 
CDF [ia and D0 

In parallel, the XenonlOO Collaboration has re- 
cently released results from a search for direct 
spin-independent dark matter scattering with 



^ Other analyses can be found in 1131 . where similar effects 
were found in the CMSSM and in gauge-mediated models. 
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100.9 live days of data using a fiducial target with 
a mass of 48 kg [21]. As we see later, this provides 
constraints on the parameter spaces of supersym- 
metric models that complement those provided 
by collider experiments Q 

In this paper we combine these new con- 
straints in updated global frequentist analyses of 
the parameter spaces of the CMSSM, NUHMl, 
VCMSSM and mSUGRA that take into account 
the results of all the searches using 2010 LHC 
data as well as the new XenonlOO constraint on 
the spin- independent scattering cross section, Cp^. 
At each point in the parameter spaces of these 
models, we construct a global likelihood function 
using previous data on electroweak precision ob- 
servables, (g — 2)^ and BR(6 — ?► sj), and applying 
the strongest of the new constraints from searches 
for multijet +]^t events, in combination with the 
constraints from H/A searches, BK{Bs — fi'^fj,~) 
and (Tp^, via the implementations described in the 
next section. 

The ATLAS and CMS searches for multijet -|- 
events provide constraints in complementary 
regions of the (mo,mi/2) planes of these mod- 
els, while the searches for heavier neutral MSSM 
Higgs bosons provide a relevant constraint in the 
(Ma, tan /3) plane of the NUHMl. The LHCb 
search for BK{Bs — >■ in combination 

with the CDF and D0 searches, affects signifi- 
cantly the likelihood function for this observable, 
with particular relevance for the NUHMl. The 
best-fit points in our new fits including all these 
2010 LHC constraints and the limit from the 
XenonlOO experiment are all close to or within 
the regions favoured by pre-LHC fits at the 68% 
CL. The spectra are somewhat heavier in the 
cases of the CMSSM, NUHMl and VCMSSM, 
whereas the best-fit mSUGRA spectrum is lit- 
tle changed. The XenonlOO upper limit on tTp^ 
has little impact on the favoured regions of the 
VCMSSM and mSUGRA, and the impact on the 

^See 1251 for discussions of the XenonlOO results in the 
context of various models including the CMSSM. Ref. 1261 
compares LHC limits and the sensitivities of astrophys- 
ical searches for supersymmetric dark matter in spe- 
cific CMSSM ('Tio,'"i/2) planes for fixed values of tan/3. 
Ref. 1271 discusses the interplay between XenonlOO and 
LHC searches in the context of a no-scale flipped SU(5) 
model. 



CMSSM and NUHMl parameter spaces is limited 
by the present experimental uncertainty in the 
hadronic scattering matrix element, that is cur- 
rently inherited primarily from the uncertainty 
in the low-energy tt-N a term, S^at. Based on 
the combination of 2010 LHC and XenonlOO con- 
straints, we present updated likelihood functions 
for sparticle masses and other observables includ- 
ing rrig, BR(i?s — > l-i-^l-i^) and ct^^. We also 
present predictions for the spin-dependent scat- 
tering cross section, (7^°, that lie considerably 
below the present experimental upper limits. Fi- 
nally, as an offshoot of our analysis, we discuss 
briefiy the potential impact of our results on fu- 
ture e'^e" colliders. 

2. Methodology 

Our analyses are performed using the 
MasterCode framework H [HHHl UHl US ■ The 

analyses have been made in a frequentist ap- 
proach, in which we construct a global likeli- 
hood function with contributions from precision 
electroweak observables, B-physics observables, 
{g — 2)^ and the astrophysical cold dark mat- 
ter density Qy^h^ as well as the limits from 
the direct LEP searches for the Higgs boson 
and sparticles and, most recently, from sparti- 
cle searches at the LHC. The model parameter 
spaces are sampled using Markov Chain Monte 
Carlo (MCMC) techniques described in our pre- 
vious papers. Our previous MCMC samplings of 
the CMSSM and NUHMl parameter spaces each 
comprised some 25,000,000 points, whereas those 
of the VCMSSM and mSUGRA include some 
30,000,000 and 17,000,000 points, respectively. 
For the purposes of this paper we have added a 
sample of some 5,000,000 CMSSM points with 
mo < 600 GeV and 250 GeV < mi/2 < 800 GeV, 
designed to improve our understanding of the 
global likelihood function at values of mi/2 that 
are somewhat larger than the previous best-fit 
values in our pre-LHC analysis of the CMSSM. 
This extra sampling had very little impact on our 
estimates of the best-fit points and 68 and 95% 
CL regions extracted from the evaluation, 
confirming the adequacy of our sampling in the 
parameter regions of interest. 
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The pre-LHC constraints are also treated simi- 
larly to our previous analyses, see Ref. [1| for the 
most up-to-date description. The numerical eval- 
uation within the MasterCode PIHIISHSI ITUII^ . 
combines SoftSUSY |29i Q, FeynHiggs [30ti33] . 
SuFla [3l[35], Superlso [Ml [37], a code for 
electroweak observables based on [38l|39] and 
MicrOMEGAs [40] (with DarkSUSY [41] as an op- 
tion not used in this paper), making extensive 
use of the SUSY Les Houches Accord [HIS]. The 
predictions we make for BR(i?s ^'^fj,~) using 
MasterCode are checked for specific fit parame- 
ters using the independent SSARD code [33]. In 
the analysis of in this paper, we link a part of 
SSARD to MasterCode to take account of hadronic 
uncertainties in dark matter scattering matrix el- 
ements, making cross-checks with MicrOMEGAs. 

The MasterCode is designed in such a way that 
the constraints from new observables can be taken 
into account and incorporated quickly and easily 
into the global likelihood function as 'afterburn- 
ers', i.e., by adding the calculated contribution 
to the likelihood function from the new observ- 
able and subsequently re-evaluating the global 
function. The new ingredients in this analysis 
coming from 2010 LHC and other searches are 
incorporated as just such 'afterburners', via the 
implementations described below. 

3. Implementations of 2010 LHC and other 
Constraints 

Studies by the LHC Collaborations have shown 
that multijet -I- constraints, with or with- 
out a single lepton, are relatively insensitive to 
tan/3 and Aq. Accordingly, we treat the AT- 
LAS and CMS constraints on such signatures as 
independent of tan/3 and Aq, and regard their 
constraints in the (mo,mi/2) plane as 'univer- 
sal' [lUinilTllllS]. At each point in this plane, 
we compare the strengths of these ATLAS and 
CMS constraints, and retain the stronger, not at- 
tempting to combine the constraints from differ- 
ent experiments. 

The constraints due to CMS and ATLAS 

^In this paper we have upgraded from the version 2.0.11 
used in earlier analyses to the new version 3.0. 13: we 
indicate below where this change affects our analysis. 



searches for events containing two or more lep- 
tons [T5] are in general less sensitive than the 
constraints due to events with jet + and at 
most one lepton, in the models considered here, 
and hence are not relevant for our evaluation of 
the global likelihood function. Moreover, these 
searches including leptons arc also more sensi- 
tive to the value of tan /3, as are searches using b 
tags [T7|. Since the reaches of the latter searches 
do not exceed those of the pure multijet -I- I^t 
searches, even at large tan/3 ^ 50, they also do 
not contribute to the global likelihood function 0. 

ATLAS jets +1^t + 0, 1 lepton analyses 

We treat the ATLAS analyses of events with 
multiple jets, zero or one lepton and^T (ATLAS 
OL, ATLASIL) [12l[15l[16] as follows. ATLAS re- 
ports the combined results of these searches as 
a 95% CL exclusion contour in the (?Tio,mi/2) 
plane for tan/3 = 3 and = 00. As seen in [T6] . 
the ATLAS OL analysis provides the dominant 
constraint on mi/2 for mo < 300 GeV. More- 
over, Fig. 17d of shows that the ATLAS OL 
search with the greatest impact on the parameter 
spaces of the CMSSM, NUHMl and VCMSSM 
is ATLAS search D (> 3 jets with leading pT > 
120 GeV, other jets with pr > 40 GeV, > 
100 GeV, A0(jet,j/T) > 0.4, rUes > 1000 GeV, 
^T/meff > 0.25). 

Two events were observed in ATLAS OL search 
D, to be compared with the number of 2.5 ± 
1.0 ± 0.2 events expected due to SM back- 
groundsQ We interpret this as a 'signal' of — 0.5± 
2.2 events, corresponding to a 95% CL upper limit 
of 3.8 events. This corresponds to the quoted 95% 
CL upper limit of 0.11 pb and the 35/pb of in- 
tegrated luminosity analyzed by ATLAS, and re- 
produces approximately the 95% CL contour for 
search D shown in Fig. 17d of (4^. This figure 

^We note in passing that LEP and Tevatron searches for 
sparticle pair-production also do not contribute signifi- 
cantly to the global likelihood function, whereas the LEP 
search for the lightest MSSM Higgs boson docs contribute 
significantly. 

"As mentioned above, this contour is not very sensitive to 
these choices: see the discussions in I12II15||T61 . 

These errors are due to the uncorrelated systematic un- 
certainty (including also the jet energy resolution and lep- 
ton efficiencies), the jet energy scale, and the luminosity, 
respectively. 
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also reports the numbers of events expected in 
ATLAS search D for points with various differ- 
ent values of {mQ,mi^2)- We calculate the corre- 
sponding numbers of effective deviations cToff from 
the observed 'signal', and construct a map of the 
deviations for intermediate values of (t?t-o, '^^^1/2) 
by interpolating between these values. At larger 
values of (^0,7711/2), where expected event num- 
bers are not provided, we scale the event numbers 
oc M~*, where M = ^Jm^'+^m^J^, following [5] 

and consistent with previous ATLAS studies. We 
then estimate the corresponding numbers of ef- 
fective deviations actf from the observed 'signal' 
using the same prescription as above, and use this 
to calculate the corresponding value of x^- 

For mo > 300 GeV, the best available ATLAS 
constraint on mi/2 comes from a combination of 
the ATLAS OL and ATLAS IL analyses. To esti- 
mate the corresponding contribution to the likeli- 
hood function at larger {1710,1711/2), we again use 
M~* scaling to estimate the expected numbers of 
events. 

We evaluate the overall ATLAS contribution 
to for each of the points in our samples of 
the CMSSM, NUHMl, VCMSSM and mSUGRA 
parameter spaces by combining these treatments 
of the ATLAS searches at small and large mo . 

CMS multijet + IpT analyses 

Following the initial ar analysis [TT] that we 
analyzed previously [5], results from an addi- 
tional CMS multijet -|-^t analysis have been re- 
leased (CMS MHT) [47] which has greater sensi- 
tivity in the (mo,mi/2) plane. The CMS MHT 
analysis also imposes stronger constraints in the 
(7710, 7711/2) plane than does the ATLAS combined 
analysis jl6| when mo > 600 GeV, so we now ana- 
lyze its results in more detail. The limit obtained 
in this search is very close to the median expected 
limit, corresponding to a difference between the 
numbers of events observed and expected from 
background that is negligible compared to the fJoft 
for the number of background events. We there- 
fore approximate the impact of this search outside 
its nominal 95% CL contour again by assuming 
that the number of effective a is simply propor- 
tional to the number of signal events expected 



at any given supersymmetric point, which we as- 
sume to be (X Af"'*, following |2], and we then 
calculate the corresponding penalty. 

Combining information of ATLAS and CMS 
analyses 

In our implementation of the combination of 
these constraints, for each supersymmetric point 
we compare the contributions to from the AT- 
LAS and CMS MHT searches calculated as de- 
scribed above, and retain just the larger of the 
two penalties, dropping the contribution from 
the lesser constraint. This procedure is conser- 
vative, but any non-trivial combination of the 
constraints would require an understanding of 
the common systematic uncertainties that is cur- 
rently unavailable, and would be justified only 
if the ATLAS and CMS collaborations provided 
additional information making possible more de- 
tailed modelling of their likelihood functions. 

We note in passing that both CMS and ATLAS 
have published limits on simplified models based 
on the above searches. These limits are not di- 
rectly applicable to the classes of supersymmetric 
models considered here since, for example, they 
consider cases in which m^ ^ rrig ^ m-^o and 
gluinos decay exclusively to qqxi, whereas in the 
models considered here other gluino decay modes 
are also important. 

LHC searches for H/A — > T+r~ 

The ATLAS and CMS Collaborations have 
also released the results of searches for heavier 
MSSM Higgs bosons H/A, produced mainly via 
bb H/A and decaying to r+r~ pairs [T91I20] . 
The stronger of these constraints is provided by 
the CMS Collaboration, which we implement as 
follows. The CMS Collaboration has provided 
model-independent limits on the H/A produc- 
tion cross section times t+t~ branching ratio 
(ctxBR) at the 68%, 95% and 99.7% CLs as 
functions of Ma [IH], corresponding to a one- 
dimensional contribution of 1, 3.84, and 9, 
respectively. For each fixed value of Ma, we 
assume that the penalty for other values of 
crxBR may be approximated by the functional 
form Ax^ oc (ctxBR)^'^*^'*^ normalized to unity 
on the 68% CL line and fitting the power p{Ma) 
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independently for each value of Af^ (typical val- 
ues are ^ 1.3). The existing CMS bounds on bb — > 
H/A — > r+r^ are expected to impact signifi- 
cantly only the NUHMl scenario, where relatively 
low values of Ma and high values of tan/3 lie 
within the region allowed by other constraints at 
the 95% CL. Therefore, we have evaluated crxBR 
for a representative grid of points in the NUHMl 
by using the the SM result for a{bb Hsm) US] 
modified by the effective NUHMl couplings ob- 
tained from FeynHiggs, which we also use to cal- 
culate the branching ratio for the decay to t+t^. 
A factor of two is included to take into account 
the production of the CP-even H and the CP- 
odd A boson, which have approximately the same 
production cross section and decay widths in the 
relevant parameter space. Ma > 150 GcV and 
large tan/3. We have then checked that crxBR 
for fixed Ma has a dependence ^ i&T? f3 in the 
parameter regions of interest. Using the value of 
a X BR calculated in this way for each point in the 
NUHMl parameter space, we then apply the 
penalty estimated as described above as an after- 
burner in our global fit. 

LHCb, CDF and D0 searches for Bg A*^A*" 

The paper by LHCb [21] provides 95% and 
90% upper limits on BR(i?s M^A*") of 56 and 
43 X 10~^, to be compared with the Standard 
Model prediction of (3.2 ±0.2) x 10'^. These lim- 
its are similar to the ones provided by CDF [35] 
and D0 [23], and a combination of the results 
from the three experiments provides a stronger 
constraint on BR(i?s ^ fi^ fi^). In order to make 
such a combination, we first performed approxi- 
mate studies, based on the signal and background 
expectations in each experiment, and compar- 
ing with the observed pattern of events, generat- 
ing toy experiments that reproduce their quoted 
90% CL upper limits. The toy LHCb experi- 
ment was constructed using the infomation shown 
in Table 3 of [gT]. The toy CDF experiment 
was based on the information given in Table H 
of |22| . combined with the invariant mass resolu- 
tion, normalization factors and averaged Neural 
Network efficiencies quoted in the text. In order 
to match exactly the observed 90% limit quoted 
by CDF, a small difference in the Neural Network 



efficiencies between the CMU-CMU and CMU- 
CMX channels [22] was introduced. Finally, the 
toy D0 experiment was based on Fig. 4 of [23], 
together with the invariant mass resolution and 
normalization factor quoted in the text. These 
toy experiments also reproduce the quoted 95% 
CL limits, giving some support to this approxi- 
mate treatment. The Tevatron results were af- 
terwards recomputed using the latest world aver- 
age fd/ fs = 3.71 ± 0.47 [SO], for consistency with 
the LHCb analysis. The results of the three ex- 
periments were combined using the CLs method, 
treating the error on fd/ fs and the branching ra- 
tio of i?+ — > J/ip{fi^ ii~)K~^ as systematic errors 
common to the three experiments. The combined 
likelihood function yields formal upper limits of 
20(24) X 10-^ at the 90(95)% CL: our global fit 
uses the full likelihood function calculated using 
the above experimental information to beyond 
the 99% CL. 

XenonlOO search for dark matter scattering 

Finally, we implement the constraint imposed 
by the direct upper limit on dark matter scatter- 
ing given by the XenonlOO experiment [53]. Its 
results are presented as a 95% CL upper limit 
on the spin- independent cross section as a func- 
tion of m^o , under assumptions for the local halo 
density and the dark matter velocity distribu- 
tion that arc described in |24] and have uncer- 
tainties that are small compared to that in the 
spin-independent scattering matrix element dis- 
cussed below [5ll[52]. The XenonlOO CoUabora- 
tion report the observation of 3 events in 100.9 
live days within a fiducial detector with a mass of 
48 kg, in a range of recoil energies where 1.8 ±0.6 
events were expected |f|. Using this information, 
we have constructed a model for the XenonlOO 
contribution to the global likelihood function 
as a function of the number of events using the 
CLs method, which is quite similar to a Gaussian 
function with mean 1.2 and standard deviation 
3.2 events. Our model for the XenonlOO likeli- 
hood function yields a 90% CL upper limit of 6.1 
events so, for any given value of to^o , we assume 

*The probability for such a Poisson background process 
to yield 3 or more events is 28%, so this observation does 
not constitute a significant signal. 
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that the 90% CL upper Hmit on a^^ quoted in [53] 
corresponds to 6.1 events, and use simple scaling 
to estimate the event numbers corresponding to 
other values of (t|^. We then use the Gaussian 
model for the XenonlOO function to estimate 
the contribution of this experiment to the global 
likelihood function for other (7^^ values. We note 
that, because of the insignificant 'excess' of 1.2 
events in the XenonlOO data, there is a contribu- 
tion Ax^ ^ 0.3 to the global likelihood function 
at small values of tTp^ [^. 

In order to translate this estimate into contri- 
butions to the global likelihood functions for var- 
ious supcrsymmetric models, wc must take ac- 
count of the uncertainty in the calculation of (Tp^ 
for fixed supcrsymmetric model parameters. The 
dominant uncertainty is that in the determination 
of the strange quark scalar density in the nucleon, 
(A''|ss| A'^), which is induced principally by the ex- 
perimental uncertainty in the 7r-nucleon a term, 
E^w = l/2(m„ -t- md){N\uu + dd\N): 



2{N\ss\N) 
{N\uu + dd\N) 



= 1 - 



N 



(1) 



where CTo = 1/2 (m-u -I- ind){N\uu + dd — 2ss\N) = 
36 ± 7 MeV [S3| is estimated from baryon octet 
mass splittings. Estimates of S^^r ranging from 
(To (corresponding to j/ = 0) up to a value as 
large as 64 ± 8 MeV have been given in the lit- 
erature [54] (and even larger values cannot be 
excluded [SS]), whereas a recent analysis based 
on lattice calculations [56] would suggest a lower 
value: E^jv ~ 40 MeV [57]. Here we span the 
plausible range by using as our default E^^v = 
50 ± 14 MeV, while also showing some results for 
E^AT = 64 ±8 MeV 0. 

The uncertainty in E^^r is quite significant for 
our analysis, since it corresponds to an uncer- 
tainty in the spin-independent cross section for 



^The predicted values of cTp^ at the post-2010-LHC best- 
fit points are all smaller than preferred by this 'excess', so 
they all receive Ax^ ~ 0.3 from the XenonlOO data, as 
seen in the Table [T] For this reason, the lower 68% CL 
limits on cr^'are essentially unchanged when the XenonlOO 
data are incorporated in the fits. 

^''The estimated uncertainty in ctq = 36 it 7 MeV is also 
included in our analysis, as are the smaller uncertainties 
associated with the quark masses. 



fixed supersymmetric model parameters of a fac- 
tor of 5 or more. We plea again for an effort 
to reduce this uncertainty by a new campaign of 
experimental measurements and/or lattice QCD 
calculations. 

4. Impacts of the LHC and XenonlOO Con- 
straints 

(mo,m]^/2) planes 

We display in Fig. [1] the (mo, '71-1/2) planes 
for the CMSSM (upper left), NUHMl (upper 
right), VCMSSM (lower left) and mSUGRA 
(lower right), driven by the ATLAS OL and CMS 
MHT constraints but also taking into account the 
other 2010 LHC constraints discussed above, as 
well as the XenonlOO constraint. In these and 
subsequent plots, we show in all panels best-fit 
points (in green), 68 and 95% CL regions (red and 
blue lines, respectively). Our pre-LHC results, 
taken from |2], are displayed as 'snowflakes' and 
dotted lines, and our post-2010-LHC/Xenonl00 
results are displayed as full stars and solid lines ["1. 

Pre-LHC, the most important lower limits on 
mi/2 in these models were indirect, being pro- 
vided by the lower limit on Mh from LEP, which 
had considerably greater impact in these mod- 
els than did the direct sparticle searches at LEP 
and the Tevatron. In each of the CMSSM, the 
NUHMl and the VCMSSM the direct 2010 LHC 
constraints push the best-fit values of toi/2 to sig- 
nificantly higher values, as well as their 68 and 
95% CL ranges whereas the effect of XenonlOO 
is not visible in this projection of the model pa- 
rameter spaces. Thus the direct 2010 LHC limits 
are constraining these models substantially more 
strongly than the LEP Higgs constraint. 

This can be seen explicitly in the panels of 
Fig. 121 which compare the effects of the LEP 
Higgs and 2010 LHC constraints on the CMSSM. 
The upper left panel shows the best-fit point, 
68% and 95% CL contours without applying ei- 
ther the LEP or the 2010 LHC constraints, and 



^^Our pre-LHC results differ slightly from those given 
in [2] as we use updated software including SoftSUSY 
3.0.13. 

^^On the other hand, the best-fit mSUGRA point is raised 
somewhat less, due to the different form of the global 
function. 
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Figure 1. The (mo,TOi/2) p/anes m the CMSSM (upper left), the NUHMl (upper right), the VCMSSM 
(lower left) and mSUGRA (lower right). In each plane, the best-fit point after incorporation of the 2010 
LHC and XenonlOO constraints is indicated by a filled green star, and the pre-LHC fit by an open star. 
The 68 and 95% CL regions are indicated by red and blue contours, respectively, the solid lines including 
the 2010 LHC and XenonlOO data, and the dotted lines showing the pre-LHC fits. 



the upper right panel shows the effect of apply- 
ing the LEP Higgs constraint but not the 2010 
LHC constraints. We note that LEP moves the 
best-fit from - 270 GeV to - 320 GeV 

while the 95% CL contour at large toq E^nd mi/2 
expands slightly, reflecting the small rise in the 
minimum of x^. The lower left panel shows the 
best-fit point, 68% and 95% CL contours apply- 
ing the 2010 LHC constraints without the LEP 
Higgs constraint. The best-fit point now moves 
to TOi/2 ^ 470 GcV, and the 95% CL contour 
moves correspondingly much further out. Finally, 
the lower right panel shows the effect of apply- 
ing the LEP Higgs constraint as well. We see 



that the best fit remains essentially unchanged at 
mi/2 ~ 470 GeV; and the 95% CL contour is lit- 
tle affected at large mg and mi/2- In summary, 
applying the LEP Higgs constraint increases mi/2 
by ^ 50 GeV in the absence of the LHC con- 
straints, and only marginally if they are applied, 
whereas the LHC constraints increase toi/2 by 
^ 200 GcV in the absence of the LEP Higgs con- 
straint, and by ^ 150 GeV if it is applied. Cor- 
respondingly, the effects of LHC on the 95% CL 
contour arc much greater than those of the LEP 
Higgs constraint. 

As seen in Fig. [l] the effects of the LHC on 
the best-fit values of mo are smaller, though 
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Figure 2. The {^0,^1/2) planes in the CMSSM with neither the LEP Higgs constraints nor the LHC 
constraints applied (upper left), with LEP hut without the LHC (upper right), without LEP hut with the 
LHC (lower left) and with hoth LEP and the LHC (lower right). In each plane, the hest-fit point is 
indicated hy a filled green star, and the 68 and 95% CL regions are indicated hy red and blue contours, 
respectively. 



there are significant increases in the CMSSM and 
VCMSSM that arc correlated with the increases 
in mif2- We note that in aU models the new best- 
fit point lies within or close to the border of the 
pre-LHC 68% CL contour, indicating that there 
is no significant tension between the LHC con- 
straints and prior indications on the scale of su- 
persymmetry breaking. Nevertheless, in all cases 
other than mSUGRA, the pre-LHC best fit points 
are now excluded at the 95% CL. Furthermore, 
the 2010 LHC constraints exclude roughly half of 
the pre-LHC 68% CL regions in the CMSSM and 
VCMSSM, and most of the pre-LHC 68% CL re- 
gion in the NUHMl. However, the LHC has yet 



to make any significant inroad into even the 95% 
CL region of mSUGRA[3. 

In Table [1] we compare the post-2010- 
LHC/XenonlOO best-fit points found in this pa- 
per with pre-LHC resuhs g] in the CMSSM, 
NUHMl, VCMSSM and mSUGRA (in the latter 
case, only the best fit in the coannihilation region 
is reported) . In addition to the minimum value of 

, the number of degrees of fi:eedom, and the fit 



^■^The raggcdness of the CL contours should be regarded 
as indicative of the uncertainties in our analysis. Recall 
also that, as already mentioned, our pre-LHC results differ 
slightly from those given in [2], as updated software was 
used, in particular SoftSUSY 3.0.13. 
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Model 


Minimum 


Probability 


mi/2 


mo 


Ao 


tan [3 


Ah (GeV) 




xVdof 




(GcV) 


(GcV) 


(GeV) 




(no LEP) 


CMSSM pre-LHC 


22.5/19 


26% 








lOtf 


108.6 


post-2010-LHC 


26.1/19 


13% 


470ii« 






99+27 


115.7 


post-Xenon (50 ± 14) 


26.2/20 


16% 


470+1^0 




-780t^^J° 


99+27 


115.7 


NUHMl pre-LHC 


20.5/17 


25% 




loot™ 


920t?^°o 


7+11 
' -2 


119.4 


post-2010-LHC 


24.1/18 


15% 




iin+SO 


O7n+1070 
■J '"-1000 


07+24 
^' -10 


117.9 


post-Xenon (50 ± 14) 


24.2/19 


19% 


530tfo° 




O7n+1070 
•-'"-'-1000 


07+24 
^' -10 


117.9 


VCMSSM pre-LHC 


22.6/20 


31% 


300^^° 


60t- 


on+50 
"^"-30 


81:1 


110.0 


post-2010-LHC 


27.9/20 


11% 




iin+UO 


^^"-190 




115.0 


post-Xenon (50 ± 14) 


28.1/21 


14% 


470^^^" 


iin+llO 


120t?°^ 




115.0 


mSUGRA pre-LHC 


29.4/19 


6.0% 




230t^^ 


430t^^" 




107.8 


post-2010-LHC 


30.2/20 


6.7% 


650t™o 


270t^° 


r;or)+130 
^■^"-130 




122.2 


post-Xenon (50 ± 14) 


30.3/21 


8.6% 


650tI°o 


2701^5° 


^■^"-130 


30tt 


122.2 



Table 1 



Comparison of the best-fit points found in the pre-LHC analysis in the CMSSM, the NUHMl, the 
VCMSSM and the coannihilation region of mSUGRA \^\^^, and our latest results incorporating the 
CMS, ATLAS, LHCh, CDF, D0 and XenonlOO constraints. We also include the minimum value of 
and the fit probability in each scenario, as well as the predictions for without imposing the LEP 
constraint. 



probability in each scenario, we include the val- 
ues of mi/2,rno, Aq and tan/3 at all the best- fit 
points, as well as the respective one-dimensional 
68% CL ranges, and the predictions for Mh if 
the LEP Higgs constraint is neglected. We note 
again that the 2010 LHC constraints are signifi- 
cantly stronger than those from previous sparti- 
cle searches and the LEP Higgs limit, resulting 
in significant increases in the best-fit values of 
^-nd smaller increases in toq in the CMSSM, 
NUHMl and VCMSSM. We note also significant 
increases in the best- fit values of tan/3 in these 
models, which are required by the {g — 2)^ con- 
straint in order to compensate for the larger val- 
ues of TOi/2 and mo. In the case of the VCMSSM, 
the scope for increasing tan /3 is restricted by 
the condition that Ao = i?o + "ZOj which is 
largely responsible for the relatively large increase 
in post-2010-LHClll The values of Aq are 
poorly constrained in all the models, and we have 

-'^''We recall that our convention [2][8] for the sign of Aq is 
opposite to that of Sof tSUSY. 



checked that there is not a strong dependence of 
the of the NUHMl on the non-universality be- 
tween the soft supersymmetry-breaking contribu- 
tions to the Higgs and sfermion masses, though 
small values of the former are somewhat pre- 
ferred. We see that the minimum values of 
have been increased by the inclusion of the 2010 
LHC data, in particular. These increases re- 
sult in some decreases in the overall probabilities, 
though insufiicient to call the models into ques- 
tion. The XenonlOO constraint causes only small 
changes in the best-fit parameters of the models 
studied, as well as small increases in the values 
and a corresponding small increase in the proba- 
bility. 

Since the constraint that most disfavours large 
supersymmetry-breaking masses is {g — 2)^, and 
since it is the interplay between this and the ad- 
vancing LHC constraints that pushes the best 
fits towards larger values of tan/3, we have in- 
vestigated the effect of dropping this constraint 
altogether. This possibility was explored previ- 
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ously using the prc-LHC data set in [7], where 
it was found that the largc-mo focus-point re- 
gion was shghtly disfavoured in the CMSSM and 
NUHMl, even when dropping the {g — 2)^ con- 
straint, by a combination of other observables in- 
cluding Mw, in particular. Now, when (g — 2)^ is 
dropped, using the 2010 LHC data set (whether 
the XenonlOO constraint is included, or not) we 
find a secondary minimum in the focus-point 
region that is disfavoured in the CMSSM by 
Ax^ ^ 1.0 , whereas this region was disfavoured 
by Ax^ 1.6 when (g — 2)^ was dropped from 
the pre-LHC data set. In the case of the NUHMl, 
we do not find a clear secondary minimum in the 
focus-point region when (17 — 2)^ is dropped post- 
20 10-LHC. 

(tan/?,TOx/2) planes 

In Fig. [3] we display the (tan/3, planes 
for the CMSSM (upper left), NUHMl (upper 
right), VCMSSM (lower left) and mSUGRA 
(lower right). We see again that the best-fit point 
and likelihood contours in mSUGRA arc only 
mildly affected by the LHC data, whereas there 
are significant increases in the best-fit values of 
tan/3 in the CMSSM, NUHMl and VCMSSM 
that are correlated with the increases in mi/2- As 
already commented in [2] , these increases may be 
understood from the interplay of the LHC and 
{g — 2)^ constraints. It is well known that for 
fixed tan/3 (g — 2)^ favours an elliptical band in 
the (mo,mi/2} plane that moves to larger mass 
values as tan/3 increases. Hence the pressure of 
the LHC towards larger values of mi/2 is recon- 
ciled with (g — 2)^ by increasing tan /3. It is ap- 
parent from the upper panels and the 68% CL 
ranges given in the Table that the constraints on 
the possible values of tan/3 in the CMSSM and 
NUHMl were quite weak pre-LHC, and are still 
not very strong. In the lower left panel, we see 
that in the VCMSSM the best-fit value of tan /3 
has increased and its range has broadened con- 
siderably post-20 lO-LHCEfl. 



^ Many early LHC analyses assumed tan /3 = 3 as a de- 
fault. It is apparent from Fig. [3] that such low values were 
disfavoured even pre-LHC, and that a more plausible de- 
fault choice post-2010-LHC would be tan /3 = 10 or more. 



(Myi,tan/3) planes 

We display in Fig. 2] the corresponding best- 
fit points and 68 and 95% CL regions in the 
(Ma, tan/3) planes for the CMSSM, NUHMl, 
VCMSSM and mSUGRA including the 2010 LHC 
and XenonlOO constraints. The LHC bb — > 
H/A — > T+T~ constraint has some impact in the 
NUHMl, where a small part of the upper left re- 
gion of the NUHMl {Ma, tan (3) plane has been 
disfavoured by this new constraint, whereas the 
previous Tevatron constraints on H/A produc- 
tion had not impacted significantly the parameter 
spaces of any of the models. 

Fig. El illustrates the effects of the CMS H/A 
constraint and the LHCb/CDF/D0 BR{Bs 
fi^fi") constraint on the (Ma, tan (3) plane in the 
NUHMl. The other LHC constraints are applied 
in all panels, but not the XenonlOO constraint. 
The left panel drops both the H/A — ^ r+r" and 
BR(i?s — > jX^ jjL^) constraints, and the right panel 
includes both constraints, and we note two prin- 
cipal effects. One is a contraction in the 68% 
CL region at lower Ma, resulting in the 68% CL 
lower limit on Ma increasing from ^150 GeV to 
^ 200 GeV, which is due to the H/A constraint. 
The other effect is some erosion of the 68% CL re- 
gion at large tan fi > 50, reducing the upper limit 
on Ma from ~ 600 GeV to - 550 GeV, which is 
due to the BR{Bs — > fi'^ fi^) constraint. However, 
we observe that the location of the best-fit point 
at {Ma, tan 13) - (400 GeV, 26) is quite insensi- 
tive to these constraints, indicating that they are 
not yet attacking the 'heartland' of the NUHMl 
parameter space. 

(j4o/mo,tan/3) planes 

Fig. [S] displays the (Ao/mg, tan/3) planes in 
the CMSSM (upper left panel), in the NUHMl 
(upper right panel), in the VCMSSM (lower left 
panel) and in mSUGRA (lower right panel). We 
see that the effect of the 2010 LHC constraints 
in the CMSSM is to push the preferred region 
towards negative values of Ao/rnQ, largely as 
a result of the push towards larger values of 
tan/3 required to reconcile the LHC data with 
{g — 2)^. The effects of the available constraints 
in this plane are weaker in the NUHMl, particu- 
larly for larger values of tan/3. In the cases of 
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Figure 3. The (tan/3,mi/2) planes in the CMSSM (upper left), the NUHMl (upper right), the VCMSSM 
(lower left) and mSUGRA (lower right). In each plane, the best-fit point after incorporation of the 2010 
LHC and XenonlOO constraints is indicated by a filled green star, and the pre-LHC fit by an open star. 
The 68 and 95% CL regions are indicated by red and blue contours, respectively, the solid lines including 
the 2010 LHC and XenonlOO data, and the dotted lines including only the pre-LHC data. 



the VCMSSM and mSUGRA, we see that the 
(Aq/too, tan/3) planes are quahtatively similar to 
those shown in [5], the main difference being a 
shift of the best-fit point to larger Ao/mg and 
tan j3. 

Cluino mass 

Fig. [7] illustrates the impacts of the 2010 LHC 
data on the likelihood functions for mg in 
the different models. The plots display the Ax^ 
contributions of the different fits relative to the 
respective best-fit points. The pre-LHC likeli- 
hood functions are shown as dotted lines, and 
the post-20 10-LHC likelihood functions as solid 



lines. In each of the CMSSM, NUHMl and 
VCMSSM, the general effect of the 2010 LHC 
data is to increase the preferred value of nig by 
^ 300 GeV beyond our pre-LHC analyses [5], 
reaching ^ 1000 — 1300 GeV, which is also some 
100 GeV beyond the results of our previous anal- 
yses using the initial CMS ar and ATLAS IL 
searches [5], whereas there is no significant effect 
on the likelihood function for nig in mSUGRA. 
Since the plots display the relative Ax^ contribu- 
tions, the differences in the overall between the 
pre- and post-2010-LHC minima of ^ 4(6) in the 
CMSSM/NUHMl (VCMSSM) are responsible for 
the differences between the pre- and post-2010- 
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Figure 4. T/ie (MA,tan/3) planes in the CMSSM (upper left), the NUHMl (upper right), the VCMSSM 
(lower left) and mSUGRA (lower right). In each plane, the best-fit point after incorporation of the 2010 
LHC constraints is indicated by a filled green star, and the pre-LHC fit by an open star. The 68 and 95% 
CL regions are indicated by red and blue contours, respectively, the solid lines including the 2010 LHC 
data, and the dotted lines including only the pre-LHC data. 



LHC likelihood functions at large nig, where the 
LHC constraints have no effect on the absolute 
values of the functions. These new normal- 
izations of are also responsible for the appear- 
ances of high-lying secondary minima at nrig 
400 GeV in the CMSSM and VCMSSM, respec- 
tively, which was previously located out of sight 
at Ax^ > 9 for the CMSSM but has now dropped 
into view. These secondary minima, like that for 
mSUGRA, are compatible with the astrophysi- 
cal cold dark matter density constraint thanks to 
rapid annihilation through a direct-channel light 
h pole. The primary minima are located in the 
fi — x5 coannihilation regions, whereas the focus- 



point regions arc strongly disfavoured in our anal- 
ysis, and not seen in any panel of Fig. [T] 

Fig. [5] displays the post-2010-LHC likelihood 
functions for BR(_Bs fi~^fi~), normalized to 
the SM prediction, where we see two principal 
effects. In the CMSSM and, to some extent, also 
in the VCMSSM (upper and lower left, respec- 
tively), values of BR{Bs pL^ exceeding the 
SM prediction are less disfavoured than in the 
pre-LHC case. This effect has a twofold origin. 
On the one hand, the LHC data disfavour a re- 
gion of parameter space where a negative inter- 
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Figure 5. The (Af^,tan/3) planes in the NUHMl with neither the CMS H/A — s> r+r^ constraint nor the 
BR(_Bs —5- jjL^ constraint applied (left) and with both the H/A and the BR(i?s — fi'^fi~) constraints 
(right). In each plane, the best- fit point is indicated by a filled green star, and the 68 and 95% CL regions 
are indicated by red and blue contours, respectively. 



ference between SM and non-SM amplitudes gave 
rise to BR(i?s — >• slightly below the SM 

prediction. On the other hand, the LHC data 
increase the minimum of significantly, and in 
this figure we show the Ax^ contribution relative 
to the respective best-fit point. Since the abso- 
lute values of a-t large BK{Bs — !> arc 
essentially unchanged by the LHC data, the dif- 
ference between the values of at the minimum 
and at large BR(i3s — > are also reduced 

by ~ 4(6) in the CMSSM/NUHMl (VCMSSM). 
In the NUHMl case we observe another impor- 
tant effect: here values of BR(_Bs — > much 
greater than the SM value (by a factor more than 
about 6) arc now more disfavoured than in our 
previous analysis. This is due to the implemen- 
tation of the LHCb, CDF and D0 constraints on 
BR{Bs — > that increases substantially the 

values at large BR(i?s Neverthe- 
less, we stress that a value of BR(i?s — > ^'^fi~) 
that is substantially larger than the SM value is 
still more likely in the NUHMl than in the other 
models (note the different horizontal scale used 
for the NUHMl). 

Light Higgs mass predictions 

In Fig.[5]the one-parameter functions for the 
lightest MSSM Higgs mass Ah in the CMSSM, 



NUHMl, VCMSSM and mSUGRA are shown. In 
this figure we do not include the direct limits from 
LEP [SHllSn] or the Tevatron, so as to illustrate 
whether there is a conflict between these lim- 
its and the predictions of supcrsymmctric mod- 
els. For each model we display the new likeli- 
hood functions corresponding to the post-2010- 
LHC data set, indicating the theoretical uncer- 
tainty in the calculation of Mh of ^ 1.5 GeV by 
red bands. We also show, as dashed lines with- 
out red bands, the central value of the pre-LHC 
results (also discarding the LEP constraint). 

One can sec that in the CMSSM, VCMSSM 
and mSUGRA the heavier preferred spectra of 
the post-2010-LHC fits result in somewhat higher 
best-fit predictions for Mh. Only in the NUHMl, 
where the minimum was very shallow pre-LHC, 
docs the best-fit value come out slightly lower. 
Now all four models predict, excluding the LEP 
constraint, best-fit values for Mh above the SM 
LEP limit of 114.4 GeV [58l[59]. One other sig- 
nificant effect of the 2010 LHC data on the one- 
parameter x^ function in the NUHMl is seen in 
the region Mh < 110 GeV. Wc recall that in the 
NUHMl the LEP constraint is weakened at low 
Mh because the Z — Z — h coupling may be re- 
duced: the 2010 LHC data help to close this loop- 
hole. Now most of the preferred Mh region in the 
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Figure 6. The (Aq/too, tan /?) planes in the CMSSM (upper left panel), in the NUHMl (upper right panel), 
in the VCMSSM (lower left panel) and in mSUGRA (lower right panel). In each plane, the best-fit point 
after incorporation of the 2010 LHC and XenonlOO constraints is indicated by a filled green star, and 
the pre-LHC fit by an open star. The 68 and 95% CL regions are indicated by red and blue contours, 
respectively, the solid lines including the 2010 LHC and XenonlOO data, and the dotted lines including 
only the pre-LHC data. 



NUHMl is indeed above ~ 114 GeV, a tendency 
that was visible already in [2] . 

Spin-independent dark matter scattering 

As a preface to discussing the importance of 
the uncertainties in the hadronic matrix ele- 
ments used in the calculation of a^^, we first 
display results that ignore these uncertainties. 
In Fig. [TOl we show our previous pre-LHC, pre- 
XenonlOO results in the (m^OjfTp^) plane assum- 
ing E^Tv = 50 MeV as dotted curves, and post- 
LHC but still pre-XenonlOO results (again assum- 
ing Sttat = 50 MeV) as dashed curves (red for 



68% CLs and blue for 95% CLs), as calculated 
using SSARD [44j . We also show the corresponding 
predictions with the higher value SttW = 64 MeV 
as duller coloured curves. The current XenonlOO 
results were not used in making these predictions, 
and we display separately the 95% CL limit on 
the cross section as a function of m^a as well as 
the sensitivity bands from [24]. We see three im- 
portant effects in these plots. One is that the 
2010 LHC results push the predicted region in 
the (m^o,tTp^) plane to higher masses, but not to 
very much lower values of cr'^^ . The second effect 
is that the new XenonlOO constraint intersects the 
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Figure 7. The likelihood functions of nig in the CMSSM (upper left), the NUHMl (upper right), the 
VCMSSM (lower left) and mSUGRA (lower right). The dashed curves are derived from the pre-LHC data 
set, and the solid curves include all the 2010 LHC constraints. In each case the value of ISrx^ relative to 
the respective best-fit point is displayed. 



regions favoured in our pre- and post-2010-LHC 
analyses of the CMSSM and NUHMl. The third 
effect is that of the value of S^tv; which changes 
the predicted range of a^^ by a factor 3. The 
combination of these two latter effects means that 
any combination of accelerator and XenonlOO re- 
sults must take careful account of the uncertainty 
in TfTrN. 

We now discuss the combination of the LHC 
and XenonlOO constraints in the {m^o,ap^) 
planes when the uncertainties in the hadronic 
matrix element SttW are included, as shown in 
Fig. [TTJ As usual, the dotted lines are pre- 
LHC and XenonlOO, the dashed lines are post- 
2010-LHC but pre-XenonlOO, and the solid lines 



incorporate also the XenonlOO constraint, with 
our default assumption S^tv 50 ± 14 MeV 0. 
In the absence of the XenonlOO constraint, the 
LHC would have allowed values of (j^^ as large as 
- 10^"*^ cm2 at the 95% CL in the CMSSM and 
NUHMl, as seen in Fig. [TOl whereas only values 
below ~ 10"^''^ cm^ would have been expected 



^^These planes cannot be compared directly to those 
in [2l|8], because here we use the SSARD code 1441 to 
evaluate cr^'. This allows a more complete treatment 
of different contributions to the scattering rates than 
does MicrOMEGAs, including important uncertainties in the 
hadron scattering matrix elements 1521 . These lead, in par- 
ticular, to larger ranges of cr^' for fixed values of m -o . We 

note in passing that MicrOMEGAs [40] uses S^jv = 55 MeV 
as a default. 
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Figure 8. The likelihood functions of BR(i3s — fi ) relative to the SM prediction in the CMSSM 
(upper left), the NUHMl (upper right), the VCMSSM (lower left) and mSUGRA (lower right). The 
dashed curves are derived from the pre-LHC data set, and the solid curves include all the 2010 LHC 
constraints. In each case the value of A^^ relative to the respective best-fit point is displayed. 



at the 95% CL in the VCMSSM and mSUGRA. 
Since XenonlOO imposes < 5 x 10~^* cm^ 
for ruj^o ^ 200 GeV, this constraint has signifi- 
cant impact in the CMSSM and NUHMl, as one 
could expect. 

In Fig. [12] we compare our predictions for cr^ 
after incorporation of the 2010 LHC data set and 
the XenonlOO constraint, for two different choices 
of S^jv = 50 ±14 MeV (our default choice, shown 
in brighter colours) and 64 ±8 MeV (a less conser- 
vative choice, shown in duller colours). As usual, 
the upper left panel shows predictions for the 
CMSSM, the upper right panel for is the NUHMl, 
the lower left panel shows the VCMSSM, and the 
lower right panel is for mSUGRA, and the 68% 



(95%) CL regions are indicated by solid red (blue) 
contours. In all models, we find that the upper 
limits on are rather independent of the value 
assumed for Sttat. However, the lower bounds 
on (jp^ are quite different for our default assump- 
tion S^AT = 50 ± 14 MeV and the comparison 
choice Yjt,m = 64 ± 8 MeV, differing by a factor 
3. This means the interpretation of future di- 
rect dark matter search constraints will be ham- 
strung by this uncertainty. 

Spin- dependent dark matter scattering 

The SSARD code also provides as an output the 
spin-dependent LSP-proton cross section, , 
and we display in Fig. [T3] the predictions for 
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Figure 9. The one-parameter likelihood functions for the lightest MSSM Higgs mass Mh in the CMSSM 
(upper left), NUHMl (top right), VCMSSM (lower left) and mSUGRA (lower right). In each panel, 
we show the functions of the post-20 10-LHC/Xenonl 00 constraints as solid lines, with a red band 
indicating the estimated theoretical uncertainty in the calculation of Mh of ^ 1.5 GeV, and the pre-LHC 
function is shown as a dashed line. 



CTp from our likelihood analysis. We see that 
the range of (t|^ is much wider in the NUHMl 
than in the other models, with both larger and 
smaller values being possible. Apart from the 
supersymmetric model parameters and the local 
galactic dark matter density, which we fix here 
to be 0.3 GeV/cm^, the principal uncertainty in 



calculating the hadronic spin-dependent 

scattering matrix element, which is dominated 
by the error in the strange axial-current matrix 
element, which we take to be (A^|s7ps|A^(s)) = 
— (0.09 ± 0.03) X s^j, where is the nucleon spin 
vector. Proportionally, the uncertainty induced 
in a?^ is far smaller than that induced in a?^ 
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Figure 10. The correlation between the spin-independent dark matter scattering cross section tTp and to^o 
prior to the inclusion of the current XenonlOO results in the CMSSM (upper left panel), in the NUHMl 
( upper right panel), in the VCMSSM (lower left panel) and in mSUGRA (lower right panel). In each panel, 
we show the 68 and 95% CL contours (red and blue, respectively), the dotted curves correspond to our pre- 
2010-LHC results, and the solid lines include the 2010 LHC results. Results assuming SjrAr = 50 McV are 
shown as brighter coloured curves and E^at = 64 MeV as duller coloured curves, in each case disregarding 
uncertainties. The green 'snowfiakes' (open stars) (filled stars) are the best-fit points in the corresponding 
models. Also shown is the 90% CL XenonlOO upper limit and its expected sensitivity band. 



by the error in S^at. As we see in Fig. [T31 the 
most stringent direct experimental upper Hmit on 
CT^D due to the COUP? Collaboration ^ (sohd 
black line) lies above 10^"^^ cm^, significantly 
higher than our predictions in any of the CMSSM, 
NUHMl, VCMSSM and mSUGRA. More strin- 
gent upper limits on (7^^ are sometimes quoted 
based on experimental upper limits on energetic 
solar neutrinos that could be generated by LSP 
annihilations inside the sun [HI]. These upper 
limits often assume that the LSPs are mainly cap- 



tured by spin-dependent scattering, which is not 
the case in general, and are in equilibrium in- 
side the sun, which is also not the case in gen- 
eral j62| . They also make simplifying assumptions 
about the annihilation final states that are not in 
general valid in the specific models studied here. 
Even with these assumptions, the upper limits 
lie above the ranges we predict in the CMSSM, 
VCMSSM and mSUGRA and barely touch the 
NUHMl range. Therefore, a direct confrontation 
of these models with data on energetic solar neu- 
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Figure 11. The correlation between m^o and the spin-independent dark matter scattering cross section 
calculated assuming a n-N scattering a term S^at = 50 ± 14 MeV in the CMSSM (upper left panel), in 
the NUHMl (upper right panel), in the VCMSSM (lower left panel) and in mSUGRA (lower right panel). 
In each panel, we show as solid (dashed) lines the 68 and 95% CL contours (red and blue, respectively) 
after (before) applying the XenonlOO \24^ constraint. The green filled (open) stars are the best- fit points 
in each model including (excluding) the XenonlOO data. Also shown are best fit and 68 and 95% CL 
contours obtained from the pre-2010-LHC data set excluding the XenonlOO result (green 'snowflake' and 
dotted lines). 



trinos still lies in the future. 

5. Summary and Discussion 

We have explored in this paper the implications 
of the 2010 LHC data for some of the simplest re- 
alizations of the MSSM, namely the CMSSM, the 
NUHMl, the VCMSSM and mSUGRA. In addi- 
tion to the most sensitive available ATLAS and 
CMS searches for jets -I- I^t, we have incorpo- 
rated the constraints imposed by searches for the 



heavy MSSM Higgs bosons H/A — > r+r~ and 
the constraints imposed by LHCb, CDF and D0 
on BR{Bs — )■ and we have also explored 

the impact of the direct XenonlOO search for dark 
matter scattering. 

We have found that the ATLAS OL and CMS 
MHT analyses shift the preferred regions in the 
(mo,mi/2) planes as compared to the situation 
based on the initial CMS ar and ATLAS IL 
searches by amounts similar to those observed 
when comparing the results incorporating the 
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Figure 12. The correlation between the spin-independent dark matter scattering cross section (7^ and 
rUj^o after including the current XenonlOO results in the CMSSM (upper left panel), in the NUHMl 
(upper right panel), in the VCMSSM (lower left panel) and in mSUGRA (lower right panel). In each 
panel, we show the 68 and 95% CL contours as solid red and blue lines, respectively. Results assuming 
^ttN = 50 ± 14 MeV are shown as brighter coloured curves, and those for E^at = 64 ± 8 MeV are shown 
as duller coloured curves. The green filled (open) stars are the best- fit points in each case. 



initial CMS ar and ATLAS IL searches with 
the pre-LHC situation. As a consequence, the 
preferred value of nig has been shifted upwards 
to 1 TeV and beyond in the CMSSM, NUHMl 
and VCMSSM. On the other hand, the picture 
in mSUGRA is not changed significantly by the 
newer ATLAS and CMS searches. 

The CMS limits on heavy Higgs production 
and our compilation of LHCb, CDF and D0 con- 
straints on BR(_Bs ^~^fj,~) have impacts on 
the parameter spaces of the NUHMl, but do not 
affect significantly the favoured regions of the 
CMSSM, VCMSSM and mSUGRA. 



The XenonlOO results have an impact on the 
model parameter spaces that would be significant 
if EjrJV were large, ~ 60 MeV. However, the cur- 
rent uncertainty in Sttw does not permit a strong 
conclusion to be drawn, and we emphasize again 
the importance of experimental and theoretical 
attempts to reduce this uncertainty. 

The adventure of the LHC search for SUSY has 
only just begun in 2010. The negative results of 
the searches to date are not in serious tension 
with the ranges of parameter spaces favoured pre- 
LHC in the models we have studied. The favoured 
regions yet to be explored offer good prospects 
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Figure 13. T/ie correlation between the spin- dependent dark matter scattering cross section ( calculated 
assuming {N\s^^s\N{s)) — —(0.09 ± 0.03) x s^, where s^ is the nucleon spin vector) and m^o in the 
CMSSM (upper left panel), in the NUHMl (upper right panel), in the VCMSSM (lower left panel) and 
in mSUGRA (lower right panel). In each panel, we show the 68 and 95% CL contours (red and blue, 
respectively) including 2010 LHC data before and after applying the XenonlOO constraint (solid and 
dashed lines, respectively). The green filled (open) stars are the best-fit points obtained with these data 
sets. Also shown are best fit and the 68 and 95% CL contours obtained from fits to the pre-2010-LHC 
data set excluding the XenonlOO result ('snowflake', dotted lines). We also show in each panel (solid 
black line) the 90% CL upper limit on provided by the COUPP Collaboration \60^ . 



for the SUSY searches during the LHC run in 
2011/12. However, it is worthwhile to consider 
whether the exclusion by the LHC of very light 
squark and gluino masses may already have mes- 
sages for future experimental studies of supersym- 
metry (if it exists). 

There is much discussion about the possible 
next large collider project to follow the LHC, with 
high-energy lepton colliders among the favourites. 
A key question is the centre-of-mass energy of 



such a collider, and indications from the LHC 
are eagerly awaited. Any definitive statement 
on the impact of LHC results must surely wait 
at least until the end of the 2011/12 LHC run, 
and will require analyses that are less model- 
specific than the results presented up to now. 
In this respect it has to be kept in mind that 
the LHC searches are mainly sensitive to the 
production of coloured particles, whereas lepton 
colliders will have a high sensitivity in particu- 
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lar for the production of colour-neutral states, 
such as sleptons, charginos and neutralinos (and 
high-precision measurements furthermore provide 
an indirect sensitivity to quantum effects of new 
states). In this sense anything inferred from the 
coloured sector on the uncoloured sector depends 
on the underlying model assumptions, and in par- 
ticular on assumptions about a possible universal- 
ity of soft SUSY-breaking at the GUT scale. 

The upward shifts for the preferred values of 
TOi/2 and, to a lesser extent, toq, that we have 
found in the CMSSM, NUHMl and VCMSSM 
upon inclusion of the 2010 LHC and XenonlOO 
constraints translate within those models into 
corresponding shifts in the production thresholds 
of supersymmetric particles at e^e~ colliders. It 
has to be noted, however, that together with this 
upward shift of the preferred mass values we ob- 
serve a significant decrease in the fit probabilities 
of those simple models, see Table [1] This indi- 
cates a slight tension in those models between the 
preference for rather light colour-neutral states 
arising in particular from {g — 2)^ and the search 
limits from the direct searches for coloured SUSY 
particles at the LHC. The mSUGRA scenario 
yields a significantly worse description of the data 
than the other considered models already for the 
pre-LHC data set, and inclusion of the 2010 LHC 
and XenonlOO constraints has only a small im- 
pact on the preferred fit values and the fit prob- 
abilities. If the upcoming LHC results lead to 
a further increase of the excluded mass regions 
for coloured superpartners, the CMSSM, NUHMl 
and VCMSSM scenarios could eventually get un- 
der pressure. Such a tension could be avoided in 
realisations of SUSY with a larger splitting be- 
tween the coloured and the colour-neutral part of 
the spectrum (for instance in GMSB-type scenar- 
ios), such that the masses of squarks and gluinos 
are in the TeV range, while sleptons, neutralinos 
and charginos can still be light. 

Alternatively, the spectrum could be com- 
pressed, decreasing the splitting between the 
coloured and colour-neutral spartieles, leading to 
softer jets from gluino and squark decays, and 
hence less stringent constraints form searches for 
jets -I- missing transverse energy at the LHC [55] . 



Noted Added 

The analysis of this paper provides a baseline 
with which 2011 LHC data can be confronted. 
While completing this work, we became aware 
of preliminary results from an analysis of events 
with > 2 jets, missing transverse energy and no 
detected leptons obtained with 165/pb of 2011 
ATLAS data [64]. In Fig. [14] we superpose on 
the (mo,mi/2) planes shown previously in Fig. [T] 
the preliminary 95% CL limits obtained using 
a PCL approach (solid black line) and a CL^ 
approach (dash-dotted black line). We see that 
in the CMSSM (upper left) and VCMSSM (lower 
left) the new preliminary PCL 95% contour runs 
very close to the best-fit point we find with the 
combined 2010 LHC data. It runs somewhat 
further away from the NUHMl best-fit point, 
and outside our 68% CL contour for mSUGRA. 
The preliminary CL^ contour runs further below 
our best-fit points, still through our 68% CL re- 
gions for the CMSSM, NUHMl and VCMSSM 
but below our 68% CL region for mSUGRA. We 
defer incorporating this result into our analysis 
until a final version is published that enables its 
contribution to the global likelihood function to 
be modelled. However, this preliminary result 
already highlights the potential of the 2011 LHC 
run to probe deeper into supersymmetric param- 
eter space. 
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Figure 14. The (mo,mi/2) planes in the CMSSM (upper left), the NUHMl (upper right), the VCMSSM 
(lower left) and mSUGRA (lower right) as shown in Fig[^ now superposed by the preliminary 95% CL 
limits obtained by the ATLAS Collaboration \64^ using a PCL approach (solid black lines) and a CLg 
approach (dash-dotted black lines). 
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